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The long-time behavior of the survival probability for unstable multilevel systems that follows 
the power-decay law is studied based on the A'^-level Friedrichs model, and is shown to depend on 
the initial population in unstable states. A special initial state maximizing the asymptote of the 
survival probability at long times is found and examined by considering the spontaneous emission 
process for the hydrogen atom interacting with the electromagnetic field. 

PACS numbers: 03.65.Db, 32.80.-t, 42.50.Md, 42.50.Vk 



One of the crucial characters of unstable systems is the 
famous exponential-decay law. Observations of the law 
were made for many quantum systems, and its theoreti- 
cal description also proved to be attributed to the poles 
on the second Riemann sheet of the complex energy plane 
However, the deviation from the exponential decay 
law was also predicted both for short times and for long 
times 1^. Indeed, despite an apparent difficulty ^3], a 
nonexponential decay law at short times was successfully 
observed Q- On the other hand, the long-time devia- 
tion has still not been detected, even though expected 
for all systems coupled with the continuum of the lower- 
bounded energy spectrum. The main reasons behind the 
matter could be ascribed to too small survival probabil- 
ity, that is, the component of the initial state remaining 
in the state at long times. 

The unstable systems are described by the Friedrichs 
model [Ml^, which enables us to investigate the time evo- 
lution involving such processes as the spontaneous emis- 
sion of photons from the atoms 0, |1| and the photode- 
tachment of electrons from the negative ions |§.lflL ll(lll]] |. 
In the former often only the first excited level is counted, 
while other higher ones are neglected, and in the latter 
the negative ion is assumed to have only one electron 
bound state. These single-level approximations (SLA) 
could be verified as long as the lowest level is quite sepa- 
rate from the higher ones. However, the multilevel treat- 
ment of the model gives us another advantage: the choice 
of coherently superposed initial-states extending over 
various levels. In fact, it can yield a variety of temporal 
behavior that is never found in the SLA [ialTHTl[Tl|. 
Such multilevel effects on temporal behavior are still not 
well studied except for Refs. [H El El El El, and 
much less examined with respect to nonexponential de- 
cay at long times. 

In the present article, we consider the long-time behav- 
ior of the survival probability S{t) by examining the A^- 
level Friedrichs model. In particular, restricting ourselves 
to the weak coupling case, we clarify how the asymptote 
of S{t) depends on the initial states. By choosing the 



initial state localized at the lowest level, we look at the 
SLA from a multilevel treatment. Then, the result in the 
iV-level model turns out to agree with that in the SLA 
in the weak coupling regime. Furthermore, among the 
various initial states, we can find a special one that max- 
imizes the asymptote of S{t) at long times. Initial states 
that eliminate the first term of the asymptotic expansion 
of S{t) are also obtained. For clarity of discussion, we 
assume all form factors to vanish at zero energy. How- 
ever, the existence of such special initial states is proved 
to be quite general and independent of other details of 
the form factors. 

The A^-level Friedrichs model describes the couplings 
between the discrete spectrum and the continuous spec- 
trum. The Hamiltonian of the model is defined by 

H^Ho + XV, (1) 

where Hq denotes the free Hamiltonian 

Hq = ujnln) {n\ / duj uj\ll!){uj\, (2) 
«=i •'0 

and AV^ the interaction Hamiltonian 

V = Y. dio [v:{Lu)\i^){n\+Vn{u;)\n){u;\], (3) 

with the coupling constant A. The eigenvalues w„ of Hq 
were supposed not to be degenerate, i.e., ujn < Wn' for 
n < n' . Both \n) and |a;) are the bound and scattering 
eigenstates of Hq, respectively, and satisfy the orthonor- 
mality condition: {n\n') — Snn', {^\^') — S{lu — lu'), and 
(n|w) = 0, where Snn' is Kronecker's delta and S{uj—u!') is 
Dirac's delta function. They also compose the complete 
orthonormal system with the resolution of identity. In 
Eq. Vn{uj) denotes the form factor characterizing the 
transition between \n) and In the latter discussion, 
we will simplify the model with the assumption that the 
form factor Vn{i^) is an analytic function in a complex 
domain including the cut (0, oo), and behaves like 



2 



where p„ and r„ are the positive constants, while qn 
and s„ are appropriate ones. The small-energy con- 
dition ensures that the integral /p (iti;u„(w)u,*/(ijj)/ti^ 
is definite. The large-energy condition ensures that 
dujVn{Lo)v*^, (w) / {z—Lo) IS definite for all complex num- 
bers z ^ [0, co). Both of the conditions are satisfied by 
several systems involving the spontaneous emission pro- 
cess of photons B B and the photodetachment process of 
electrons 0, |^ liol lll| . Note that this small-energy con- 
dition excludes the photoionization processes associated 
with the Coulomb interaction [T^ : however, the formu- 
lation developed below could be applied to those cases. 

The initial state of our interest is an arbitrary su- 
perposition of the unstable states 



N 

E 



Cn\n), 



(5) 



where c„'s are complex numbers satisfying the normal- 
ization condition = 1- Then, the survival 
probability S{t) of the initial state |?/'}, that is, the prob- 
ability of finding the initial state in the state at a later 
time i, is defined by S{t) = \A{t)\'^. The A{t) denotes 



the survival amplitude of i.e., A{t) 



-itH 



In general, the Hamiltonian Q has the possibility of pos- 
sessing not only the scattering eigenstates but also 
the bound eigenstates We shall here restrict our- 
selves to studying the decaying part of A{t), and merely 
call it the survival amplitude with the same symbol as 



A{t) 



(6) 



In order to estimate the long time behavior of A{t), let us 
evaluate the scattering eigenstates l^/'^j'^^) by solving the 
Lippmann-Schwinger equation, i.e., = + [uj ± 

iO - Hoy^XV\il^\^'^). In the case of our Hamiltonian, 
this equation can be solved in the form, iV'i^''} = |^) + 



E 



N 



1") + /o°° '^^ 



the integrand of A{t) reads, 



from which 



N 



n=l 

The F^\ijj) is determined by an algebraic equation 

N 

G„^(c. ± iO)Fi^\u:) = -Xvn{io), (8) 



where 



= {i^n - z)Snn' + A^S„„/(z), 



(9) 



which is the (n, r7,')-th component of the N x N matrix 
G^^{z), and s„„'(z) is defined by 



,{z) 



duj 



,u„(cj')<'(w') 



(10) 



for all z 



re 



tip 



{r > 0, < (fi < 2tt). Under the large- 



energy condition of Eq. Q), Snn'{z) is guaranteed to 
be analytic in the whole complex plane except the cut 
[0, oo). For the later convenience, G~^{z) is defined as 
an inverse of G{z), where G{z) is assumed to be regular. 
Note that G{z) is nothing more than the reduced (or 
partial) resolvent Gnn'{z) = {nHH — z)~^\n'). One can 
confirm this fact by following the discussion in section 
3.2 of Ref. 6]. Since the behavior of A{t) at long times 

is characterized by that of Fn^\u!) in Eq. ITJ at small 
energies, we need to estimate the small-energy behavior 
of G{z). Note that under the condition Q we have 

[/„„' (w) =F i7ru„ iuj)v*, (w)] 
+ A2/„„,(0) + o(1), (11) 



as ^ +0, where s„„'(a;±iO) = (a;) =Fz7rv„(cj)i;*, (a;) 
and /„„'(w) 



P r°° rf^^ ^-. (^') ^j^g^^ p ^^^^^ 



the principle value of the integral. The existence of 
-^nn'(O) niay be just guaranteed by the small-energy con- 
dition of Eq. Q). Supposing that G„„' is of the form 



G„„'(tJ ± iO) = gnn' + o(l), 
as cij ^ +0, one obtains that 

N 



(12) 



5nn' ~ ^nniG^j^, 
m—1 
N 

= ^ gnm[l^7nSmn' + A^/„m'(0)] + o(l), (13) 



m—1 

which leads to 

1 

Qnn' — 



N 



^nn' A ^ ^ Qnmlvan' (0) 



(14) 



We solve this equation by assuming that .g„„' can be ex- 
panded for small A as 



(15) 



By substituting Eq. ifT^ into (fn|l . it follows that 

9nl ^ Snn'/uJn', 9^1 = - Inn' (0) / UJ^UJu' , (16) 

and for j > 1 



■^7171 



1 ^ 

-T.9nrn'^^mn'{0), 



(17) 



m—1 



where we have assumed that all LOn does not vanish. Note 
that g'^^, and g^^, derived here accord with at least those 
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for solvable cases, where G{z) is explicitly obtained 0, 
IT^ . We can then obtain 



(18) 



with 



where 



N 



SlL^X^J2 +0(A^), (19) 



n' = l 



Qn (Pn = P) 
iPnJ^p) ' 



(20) 



where p = niin{p„ } . With use of the g„ instead of qn , we 
extracted only the dominant part of Fn^\uj) at small lo. 

The long time behavior of A{t) can be simply obtained 
by applying to Eg. © the asymptotic method for the 
Fourier integral |l9l |. As mentioned before, the long time 
behavior is determined by the small-energy behavior of 
its integrand. By inserting Eq. H18(l into the inte- 
grand of A{t) turns out to behave asymptotically 



N 



^ ^ fn'^n 



,2P 



(21) 



as w —> -f 0. Applying the asymptotic formula for Fourier 
integrals, we obtain from Eq. (|21|l the asymptotic behav- 
ior of Eq. © reading, 



A{t) = A' 



2r(2p- 



1) 



(ii)2P 



+ 1 



N 

E 

n=l 



+ o{t-^P-^), (22) 



as t ^ oo, where i'^+'^-P = e*('i+i-p)^/2^ and T{z + 1) 

(ixa;^e~^. We can clearly perceive A(t) ^ t^^P~^, the 
power decay law. 

Using the above result, let us first consider the higher- 
level effects on the long-time behavior that starts from 
the localized initial state at the lowest level. This study 
is directed to an examination of the SLA. For such an 
initial state, i.e., c„ = 5„i, Eq- H22|) becomes 



A{t) = X 



2r(2p 



2p+l 



0{\^)]+o{t 



(23) 



where we supposed that qi ^ 0. Since there are no factors 
related to the higher levels in Eq. (|23|) , it implies that the 
long-time asymptotic behavior of A{t) could agree with 
that in the SLA for a sufhciently small A. 

On the other hand, we can find a special superposition 
of unstable states \n) that maximizes the asymptote of 
A{t) at long times. It is worth noting that its dependence 
on the initial states only appears in Eq. 122|) through the 
factor ^n=i fn'^n, which can be rewritten by an inner 
product as 



JV 



(24) 



where we have introduced an auxiliary vector defined by 



N 



(25) 



Thus, resorting to the Schwarz inequality, we see that 
the maximum of the factor (|24|l is just attained by if and 
only if IV') oc |x), i.e.. 



Cfn 



Ixll, 



(26) 



where c is an arbitrary complex number with |c| = 1. 
Therefore, preparing the initial state [ip) according to 
the above weights (|26|l . we can maximize the asymptote 
of A{t) at long times. Substituting Eq. into Eq. 

one obtains that 



Ait) 



,r(2p + i) 



i^t) 



2p+l 



;^2r(2p 



1)^ 



2p+l 



(27) 
(28) 



It should be remarked that the initial state extended over 
unstable states |n) has the possibility of increasing the 
intensity of A{t) more than a localized one would. This 
possibility may be interpreted as follows. Let us consider 
the spontaneous emission process for an atom interact- 
ing with the electromagnetic field, where \n) is identified 
with the (n-|- l)-th excited state of the atom with the vac- 
uum state of the field and luj) is the ground state with 
the onc-photon state. In this process, an initially excited 
atom makes a transition to the ground state with emit- 
ting a photon, while the atom that fell into the ground 
state can be reexcited by absorbing a photon. In the 
latter process, there are various candidates for the ex- 
cited state. Repopulation of each excited level can make 
the intensity of A{t) grow, providing that the initial state 
possesses those excited levels. However, if the initial state 
only consists of a specific excited state, the other excited 
states composing the state at a later time t are discarded 
without any contribution to A{t) This is the rea- 

son why the decay of the A{t) for extended states can be 
relaxed more than that for localized states. 

Note that the above argument also suggests the pos- 
sibility of finding another kinds of initial states that are 
coherently superposed to eliminate the factor H24|l . This 
is indeed achieved by the initial states that are orthogo- 
nal to Ix), 



(xl^) = 0. 



(29) 



In this case, the first term in the rhs of Eq. (|22|) becomes 
zero. This fact means that A{t) for such an orthogonal 
state asymptotically decays faster than t~'^P~^. 

The maximizing initial state seems to be desirable for 
an experimental verification of the power-decay law. Let 
us now discuss the value of |^(i) p for such an initial state 
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TABLE L The level-number dependence of "^^^^ |gn/i^np, 
the decay time tiv of the (A*' + l)p-state, and the transition 
time tep from the exponential to the power decay law. 

Number of levels N 1 10 50 

|iiL|^v-iv I 9^1^ Ym L28 L29 

I qi I I uj^ I 

tN (s) 1.60 X 10"*^ 3.18 X 10"^ 3.18 x 10"^ 

tep (s) 2.00 X lO"'^ 4.23 x lO"'^ 4.59 x 10"^ 



On the other hand, the factor J2n=i \Qn/^n\'^ is essen- 
tially unchanged with N, whereas tep rapidly increases 
with N (see Tabled). These facts and Eq. imply 
that A{tep) rather decreases as N increases ,^23). Hence, 
we should unfortunately conclude that the maximizing 
initial state does not provide any help for an observation 
of the power decay law for the spontaneous emission from 
a hydrogen atom. 



at long times. In particular, we shall evaluate this value 
at the time tep of the transition from the exponential to 
the power decay law. We have to know the exact values 
of both pn and <?„ for many n's; however, this require- 
ment is satisfied by considering the spontaneous emission 
process for the hydrogen atom interacting with the elec- 
tromagnetic field (see also Refs. This time, |n) 
is interpreted as the (n + l)p-state of the atom with the 
vacuum state of the field, and |w) as the Is-state with 
the one^hoton state. It then follows that p„ = 1/2 for 
all n |21| . and \qn\ is determined through the relation 
7„ = 27rA2|?;„(w„)|2 + ©(A^) ~ 27rA2|g„|2|w„|, where the 
last estimation is confirmed in the dipole approximation. 
7„ is the decay rate of the {n + l)p-state, which is esti- 
mated as 7„ ~ 8.0 X 10^ X 2^(71 -hl)n2"/9(n-h2)2"+4 
I23. From these facts it follows that 



gn 



8.0 X lO'^ X 6{n+iyn'^ 
7rfi3(n-H2)2"+4[(n-f 1)2- 



(30) 



1I6 



that \qn\ 



LUr, 



where w„ = |f][l - (n + 1)"^] with n = 1.55 x 10^ 
s~^, and we also choose A = 6.43 x 10^^. Thus, we see 
n^^ for a large n. In Table the 

numerical values of X]n=i kn/t^nP (— llxlP), the decay 
time tM {— I/ttv) of the (TV -I- l)p-state, and the time 
tep are listed for the three cases of the level numbers 
A'' = 1, 10, and 50. Here, we define tep as the maximum 
time which equates the square modulus of the asymptote 
II28II to that of the following A{t) at intermediate times 
[11111, 



In summary, we have considered the long-time behav- 
ior of the unstable multilevel systems and estimated the 
asymptotic behavior of the survival amplitude A{t) for 
an arbitrary initial state in the long-time region where 
A{t) obeys a power decay law. We have then found two 
special initial states. One of them asymptotically max- 
imizes A{t) at long times, and the other eliminates the 
first term of the asymptotic expansion of A{t) . The latter 
fact may imply that the exponent of the power decay of 
A(t) is determined by not only the small-energy behav- 
ior of the form factors but also the initial population in 
unstable states. Such relations between the initial states 
and the power decay law were studied with respect to 
the long-time behavior of wave packets, both for the free- 
particle system ,2.4] and for finite-range potential systems 
[25j- In the case of the experimental verification of the 
power decay laws, the existence of the maximizing ini- 
tial states seems preferable. This expectation is probably 
misplaced for the spontaneous emission process of a hy- 
drogen atom, however, a possibility still could remain for 
the systems allowing the photodetachment or the pho- 
toionization process. We then should require of them the 
property that both 7„ and \qnl^n\ do not decrease as n 
increases. More important states for this aim are those 
states which maximize A[t) at the transition time from 
the exponential to the power decay law. The relation be- 
tween such a maximizing state and the discussed one is 
still unclear. It will be addressed in a future issue. 



N 



A{t) ^ 



(31) 



where c„ is chosen as Eq. (|26|) . It is worth noting that 
when t 3> tfq, \A(f)\^ can approximate |cjv|'*e~*'^" be- 
cause the decay time t„ lengthens with n in this case. 
We see from Tabled that tep is much longer than tjv, so 
that tep is roughly estimated as the root of the equation. 



i4 



N 

E 



(32) 
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